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ABSTRACT 

The vacuum cosmological model on the manifold R x Mi x . . . x M„ describing the 
evolution of n Einstein spaces of non-zero curvatures is considered. For n = 2 the Einstein 
equations are reduced to the Abel (ordinary differential) equation and solved, when (A^^i =dim 
Ml, N2 = dimM2) = (6, 3), (5, 5), (8, 2). The Kasner-like behaviour of the solutions near the 
singularity tg —>■ +0 is considered [tg is synchronous time). The exceptional ("Milne-type") 
solutions are obtained for arbitrary n. For n = 2 these solutions are attractors for other 
ones, when ts — >■ -|-oo. For dim M = 10, 11 and 3 < n < 5 certain two-parametric families of 
solutions are obtained from n = 2 ones using "curvature-splitting" trick. In the case n = 2, 
{Ni, N2) = (6, 3) a family of non-singular solutions with the topology R^ x M2 is found. 
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1 Introduction 



Usually within multidimensional cosmology (see, for instance, ||l|-|T^ and refs. therein) the 
space-time is considered manifold 

M = Rx Ml X . . . X Mn, (1.1) 

where R is the time axis, one of the manifolds Mi, . . . , M„ is interpreted as 3-dimensional 
external space and the others as so-called internal spaces. In most models Mi,...,M„ 
are considered as spaces of constant Riemann curvature, or Einstein spaces. From the 
very beginning it is supposed that the usual three spatial dimensions and extra spatial 
dimensions are on the same footing and this assumption is followed from believing (of some 
physicists) that the early Universe is multidimensional. The separation of extra dimensions 
from the usual one is attained in multidimensional cosmology by the so-called dynamical 
compactification of internal spaces to unobservable sizes when their scale factors decrease to 
the lengths of the Planck order (10~^^cm) during the Universe evolution. 

The study of dynamical properties of the multidimensional Universe such as, for instance, 
compactification of internal spaces or expanding of the external space in detail demands exact 
solutions. However, in all known integrable multidimensional cosmological models the chain 
of the constant curvature spaces Mi, . . . , M„ may contain at most one space of non-zero 



curvature (see, for example, [|T0|)- As far as we know, there are no integrability conditions 
or explicit integration methods for the models describing evolution of two or more Einstein 
spaces with non-zero Ricci tensors. The aim of this paper is to integrate the vacuum model 
on the manifold R x Mi x M2, where the Einstein spaces Mi and M2 have non-zero Ricci 
tensors. 

This paper is organized as follows. In Sect. 2 we describe the model and get the equations 
of motion. We show, that their integrability by quadratures is reduced to the integrability of 
some ordinary differential equation. The latter appeared to be Abel's equation. In Sect. 3, 
we show, that for the models with (dimMi,dimM2) = (6,3), (5,5), (8,2) the Abel equation 
belongs to the integrable class described by Zaitsev and Polyanin |2T[ within discrete- 
group analysis methods. We integrate the (vacuum) Einstein equations for these three models 
(see formula (3.38)) for non-exceptional case and present the trajectories of motion on the 
scale factors configuration plane. In Sect. 4 the exceptional "Milne-type" solutions (see (4.1) 
and generalization for arbitrary n in (4.2)) are presented. The Kasner-like behaviour near 
the singularity (for ts +0) is investigated in Sect. 5. In Sect. 6 certain two-parametric 
families of solutions with n curvatures for 3 < n < 5 are presented. Some non-singular 
solutions with topology R^ x M2 are considered in Sect. 7. 



2 The model 

At first, following our previous papers 0-[|ll] "we consider here the general formalism for de- 
scription of multidimensional vacuum cosmological models. It is supposed, that D-dimensional 
space-time manifold M is defined by relation (1.1). The manifold M is equipped with the 



1 



metric 



g = - exp[2'y {t)]dt dt + ^ exp[2x*(t)]5( 



(2.1) 



1=1 



where ^{t) is an arbitrary function determining the time t. It is supposed that the manifold 
Mi (see (1.1)) for z = 1, . . . , ra is the Einstein space of dimension Ni with the metric g^'^\ i.e. 



l,...,iV„ 



(2.2) 



where Aj is constant. (For the manifold Mj of constant Riemann curvature the constant 
Aj reads: Aj = Ki{Ni — 1)). 

The non-zero components of the Ricci tensor for the metric (2.1) are the following [T^ 



R', = exp[-2^]J2Ni[x' 



i=l 



exp[-27] 



(2.3) 
(2.4) 



where indices and for i = 1, . . . ,n run over from 1 + Nj to Z]}=i Nj. 

After the gauge fixing 7 = F{x'^) vacuum Einstein equations = for A, B = 1, . . . , D — 
1 and equation Rq — 6qR/2 = may be presented as the Lagrange-Euler equations and zero- 
energy constraint correspondingly obtained from some Lagrangian L(x*,x*) = Li?(x*,x*). 
For the so-called harmonic time gauge 



7 



7o = E^« 



i=l 



this Lagrangian has the following form |^ 

1 



where 



L(x\ x') = -J2 G^Ji'x^ + - E ^^^^ exp[-2x* + 270] 
i,j=i ^ i=i 



Gij = NAj - NiNj 



(2.5) 



(2.6) 



(2.7) 



are components of the minisuperspace metric. The minisuperspace metric has the pseudo- 
Euchdean signature (—,+,...,+) P, |^. (This is a well-known property of the Hilbert- 
Einstein action P^ .) 

If all spaces Mi,...,M„ are Ricci-flat, i.e. Aj = for i = l,...,n, or only one of 
them has non-zero Ricci tensor, then equations of motion following from the Lagrangian 
(2.6) are integrable and exact solutions are obtained in |]10|. But nearly nothing is known 
about integrability of the models for two or more Einstein spaces with non-zero Ricci tensor. 
Obtaining exact solutions for such models is connected with a general problem of integrability 
of the so-called pseudo- Euclidean Toda-like systems |jlO|, |l^]. These systems are described 
by the Lagrangian of the following form 



-| n—l m 

^ = ^ E Va.z'^^' - E A^'^ exp[E <z% 

^ a,fe=0 s=l a=0 



n-1 



(2i 



2 



where all u^, A^^'^ are real constants and {i]ab) =diag(— 1, +1, . . . , +1). The Lagrangian (2.6) 
may be reduced to the Lagrangian (2.8) by some linear transformation of the configuration 
variables x^, . . . ,x"'. It should be noted, that the Lagrangian (2.8) describes both vacuum 
multidimensional cosmological models or models with multicomponent perfect fluid source 



12, 14]. The well-known 4-dimensional Bianchi-IX model (the so-called "mixmaster model") 
is also described by the Lagrangian (2.8) with n = 3, m = 6 and certain u^, A*^*' p3|] . 

There are extensive literature devoted to the integrability of the Euclidean Toda-like 
systems {{rjab) =diag(-|-l, . . . , +1) in (2.8)) and methods of their explicit integration (see, for 
instance, p5]); but, as far as we know, the integrability of the pseudo- Euclidean systems is 



not well studied yet. In our previous paper [|14| we singled out two integrable by quadratures 



classes of these systems and developed their integration. The first class contains the pseudo- 
Euclidean systems trivially reducible to the Euclidean systems. It is easy to see that under 
the conditions: Uq = for s = 1, . . . , m it follows from (2.8) that i° =const. Then the La- 
grangian (2.8) effectively "loses" negative term in the kinetic energy, so for z^, . . . , z^~^ the 
Euclidean system arises. The second class appears provided vectors u^, G -R", where 

= {uq, . . . ,u^_i), are coUinear or form the orthogonal set with respect to the symmetri- 
cal bilinear form on i?" with matrix {r]ab) (see for details [T^). But the models describing 
evolution of two or more Einstein spaces with non-zero Ricci tensors are out of these classes 
and until now we had neither any integrability conditions nor explicit integration methods 
for them. In this paper we show that at least the models with n = 2, i.e. multidimen- 
sional vacuum cosmological models with metric (2.1) on the manifold (1.1), where Mi and 
M2 are Einstein spaces of certain dimensions with non-zero Ricci tensor, are integrable by 
quadratures, and develop their integration procedure. 

So, let us consider the Lagrangian (2.6) for n = 2 and Ai, A2 7^ 0. After the following 
coordinate transformation diagonalizing the minisuperspace metric: 



X = x^ — x'^ + - In 



X2N2 



y = ap (a^iX^ + N2x'^) , (2.10) 



where a and P are defined by the relations 



a^J^^l±^, P^,^, (2.11) 
V N, + N2 ' ^ V N1N2 ' ^ ' 

we obtain the following expression for the Lagrangian (2.6) in the new variables x and y 



Li^x'-f)-Vix,y), (2.12) 



L -- 

2(3' 

where the potential V{x, y) has the form 

V{x,y) = K,exp(2«/?-^?/)[sgn(Ai) exp(2/5ix) + sgn(A2) exp(2/32a;)]. (2.13) 
Here we denoted 

V^ = --\X^Nif^\X2N2\~^\ (3i = — — , p2 = ^^ — ■ (2.14) 

2' ^ ^' ' ^ ^' ' Ni + N2 N^ + N2 ^ ' 
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We note that — (^i = 1- The corresponding zero-energy constraint reads 



V{x,y) = 0. 



The Lagrangian (2.12) leads to the following equations of motion 

X = -2/?Voexp(2a/?-^?/)[/5isgn(Ai)exp(2/5ix) +/52Sgn(A2)exp(2/?2a;) 
y = 2af3V{x,y). 



(2.15) 



(2.16) 
(2.17) 



It is easy to see, that the set of equations (2.16), (2.17) does not admit solutions with 
X = y = Wt. Indeed, such solutions are possible only if A1A2 < and /5i = ^2- But the 
latter condition is impossible in the model under consideration (see (2.14)). So, our model 
has no static solutions for any dimensions of Mi and M2. 

To integrate the equations of motion (2.16), (2.17) under the zero-energy constraint (2.15) 
we consider at first the following procedure proposed in [H . Using the relation 



y 



fy 
dx"^ 



y-xy 
xP' 



equations of motion (2. 16), (2. 17) and the zero-energy constraint (2.15) in the form 

213^ 



X 



vix,y) _ ^^2 vi^.y) 



(y/xy-l " (y')'-l' 
we obtain the following ordinary differential equation of the second order 



y 



where 



fix) 



i\2 



1 



l(A+/52 + /(x))y'+^ 



tanh(a;), if A1A2 > 0, 
coth(x), if A1A2 < 0. 



f2.18) 



(2.19) 



(2.20) 



(2.21) 
(2.22) 



After this procedure is done the following solutions may be lost: x = ±y Vt and x = Vt. 
It is easy to see, that under the zero-energy constraint the equations of motion (2.16), (2.17) 
have no solutions with x = ±y Vt. The solutions with i; = Vt appear only in the case Ai < 
and A2 < 0. So, by integration of the equation (2.20) we obtain all possible trajectories on 
the configuration plane {x,y) except the trajectory x = | ln[— /5i//52] for the model with 
Ai < and A2 < 0. If the trajectories of motion are known we may get the law of motion for 
each trajectory by integration of the equation (2.19). Thus, the problem of the integrability 
by quadratures of the equations (2.16), (2.17) under the constraint (2.15) is reduced to the 
problem of the integrability of the equation (2.20). In the paper p the equation (2.20) was 
studied by qualitative methods, but here we consider its explicit integration. 

For the considered model we have from (2.11) and (2.14) 



P1+P2 



N1-N2 
N1 + N2'' 



a 



R 

N1 + N2'' 



(2.23) 
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where 

R = R{Ni, N2) = ^NiN2iNi + N2-l). (2.24) 
We also rewrite (2.21), (2.22) as follows 

. 5^pM^. (2,25) 

exp(2x) + e 

where here and below 

e = 8182, Si = sgn(Ai), (2.26) 

2 = 1,2. 

It should be noted, that the right side of the equation (2.20) does not contain the unknown 
function y due to the factorization of the potential V{x,y). So, for the function z = y' we 
have the first order equation of the form 

z' = fo{x) + h{x)z + f2{x)z^ + Ux)z\ (2.27) 



An equation of the form (2.27) is well-known as Abel's equation (see, for instance pO| , pT|). 
There are no methods to integrate Abel's equation with arbitrary functions fo{x), fi{x), 
f2{x) and fsix). 



3 Exact solutions 



If any solution ^0(2^) to the Abel equation (2.27) is known, then by the transformation |21| 
E(x) 



u = . 

Z- Zq' 

it may be written as 
where 



E{x) = exp 



{3Ux)z^o+^Ux)^o + fi{x))dx 



uu + Flix) + F2(x)u = 0, 



Fi{x) = h{x){E{x))\ F2{x) = {3h{x)zo + ^(x)) E{x) 
For the new variable 



q = — J F2{x)dx 



the Abel equation has the so-called canonical form |^ 



- M = F{q), 

aq 



where we denoted 



F{q) 



Fijx) 
F2(x) 



(3.1) 

(3.2) 
(3.3) 
(3.4) 

(3.5) 
(3.6) 



The Abel equation in the form (3.5) was recently studied by a discrete-group analysis meth- 
ods and a number of new integrable equations were described pO| , p[] . In principle, the Abel 



5 



equation (2.20) may be written in the canonical form as its solutions y' = ±1 are known, 
but some difficulty to obtain the function F{q) in an explicit form for arbitrary parameters 
a, (3 and (3i arises. It is not hard to verify, that if one of the two following conditions holds 

-a = /3(2/5i-/32), a = /3(2/52-/3i), (3.7) 

then the function F{q) for the equation (2.20) may be obtained in an explicit form by solving 
some algebraic second order equation. Using (2.11) and (2.14) we get corresponding relations 
for the dimensions of Mi and M2. 



iVi 



4Ar, 



(31 



N2-I Ni-1 
It can be easily seen that the ffist condition holds in the following three cases: i) A^^i = 6 and 
N2 = 3; ii) A^^i = 8 and N2 = 2, iii) A^^i = 5 and N2 = 5. The second condition corresponds 
to the inverse numbering of the spaces Mi and M2. 

Under one of the conditions (3.8) the Abel equation (2.20) in the canonical form may be 
reduced to the Emden- Fowler equation of the form: (pY/dX"^ =consty™°^* [^, which may 
be easily integrated. 



3.1 (l+6+3)-model 

First, let us consider only the model with A^^i = 6 and A^2 = 3. This (1 + 6 + 3)-model is 
of most interest, because the 3-dimensional space M2 may be interpreted as our (external) 
space. Omitting technical details we give at once the final Emden-Fowler equation for this 
model. One may check that the Abel equation (2.20) for A^^i = 6 and A'"2 = 3 by the 



transformation 21 



1 + eexp[2x] 



1 + 



?,XdY 
2Y dX' 

"^Y dX 



1 I X dY 
^ ~^ Y dX 



a_ ( dY\ 2 
Y \dx) 



1 



(3.9) 
(3.10) 



where a 7^ is arbitrary constant, may be reduced to the following Emden-Fowler equation 



dx^ 



aY- 



(3.11) 



which is well-known in classical mechanics as describing the one-dimensional motion of a 
charged particle in the field of the (attractive or repulsive) Coulomb center. We remind that 
e = sgn(AiA2). So, by integrating (3.11) and using the transformation (3.9), (3.10) we obtain 
the general solution to the equation (2.20) for A^i = 6 and A^2 = 3. In a parametrical form 
the result looks as follows 



e exp[2x] = $(r, -6, d) = -6 (r^ - 6) [1 + r {g{r, 
, (1 + 35r2)exp[2x] +e 



-5) + Ci)] + \., 



y 



35r2)exp[2x] +e' 



X 



In 



r 



6 



Co 



(3.12) 
(3.13) 
(3.14) 



6 



where 



r 



r + 1 
= arctaii(T), 



5 = +1, 



(3.15) 
(3.16) 



By T we denote the parameter (time), Ci and C2 are arbitrary constants. Formulas (3.12), 
(3.13) together with y' — ±1 present the general solution to the Abel equation (2.20) for 
A'"! = 6 and N2 = 3. This may be also verified by a straightforward calculation using the 
following relation 

5 



The function (3.12) is depicted in Figs. 1, 2 for different values of 5 and Ci. The limits 
$(r) — > 00 and $(r) ^0 (for non-exeptional Ci) correspond to Kasner-like behaviour near 
the singularity (see Sec. 5 below). 
Figs. 1, 2 

Using relations (3.12), (3.17), zero-energy constraint (2.15) and the expression for the 
potential (2.13) we obtain the following relation between harmonic time and r- variable 



where 



1, 



gl 11 I 2\ 



Thus, from (3.18) we obtain 

5 = -£2 = -sgn(A2). 
The transformation inverse to (2.9), (2.10) for A^^ = 6, A^2 = 3 looks as follows 



— - i 2x + V — In 
6 I ^ 



A, 



X 



1 
6 



2Ai 

4x + y + 2ln 



A2 



2Ai 



(3.18) 

(3.19) 
(3.20) 

(3.21) 
(3.22) 



Using (2.5), (3.12), (3.14), (3.18)-(3.22) we get the following relation for the metric (2.1) 
(n = 2) 

9 = c (/i/2)-^ {-2fr'dT ^dT + /2I Ail^W + I A2|5^')}. (3.23) 
Here c 7^ is an arbitrary constant and 



/i 



fi{r,e2) 



h = /2(r,ei,e2,Ci) = 2£i£2$(r,£2,C'l) 

= -3£i (r^ + £2) [1 + r (y(r, £2) + Ci)] + e^e2 > 0. 



(3.24) 



(3.25) 



We recall that £j = sgn(Ai), i — 1,2. It should be noted that although originally c > 
since c = 2"^/^^"^ exp(C2/3), the negative c also gives us the solution to vacuum Einstein 
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equations. The trajectories on the plane of scale factors > 0,i — 1,2, corresponding to 
the solution (3.23), where 

aj = c(/i/2)"^/2|Ai|, 4 = c(/i/2)"^ IA2I, 

are depicted in Figs. 3-5 (we put c = |Ai| = IA2I = 1). The points Tj = Tj(Ci), i — 1,2,3, 
are zeros of f2- 
Figs. 3 - 5 



3.2 The (1+8+2)- and (l+5+5)-models 

Some modification of the Anzatz (3.12)-(3.14) may be used for obtaining the solutions to 
the Abel equation for other two cases (A^i, A^2) = (8, 2), (5, 5). The solutions read 



£exp[2x] = $fe(T, -5,Ci) 

, _ (l + MT^)exp[2a;] +e 
^ ~ (1 -MT2)exp[2x] +£' 



y — X In 

^ 2 



r'-S 



+ C2, 



where 



k = {N, + N2)/N2 = 3, 2, 5 for (iVi, N2) = (6, 3), (5, 5), (8, 2), 
respectively and functions = $fe(T, —S,Ci) satisfy to the equations 

S 



(1 - kSr^) $fc + 1 



t(t2 - 5) 

The solutions to (3.30) are the following: $3 = $ (see (3.12)) and 



$2(t, -5, Ci) = 1 + 25 (r" -5)+ Ci5t 
1 5 



(3.26) 
(3.27) 

(3.28) 
(3.29) 

(3.30) 



$5(r, -5,Ci) = ---iT 



(r^ - 6) [St^ -2S + 3t (r^ - S) (^(r, -6) + Ci 



(3.31) 
(3.32) 



Using relations (3.26), (3.30) zero-energy constraint (2.15) and the expression for poten- 
tial (2.13) we obtain the following relation between harmonic time and r-variable 



4x + 2 ( 1 + - ) C2 



where here 



iVi 



I A^2A2| * • 



(3.33) 



(3.34) 



^From (3.33) we obtain the relation (3.20). 
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The transformation inverse to (2.9), (2.10) is the following 



No 



X In 

2 



\2N2 



AliVi 



+ 



N, + N2 



X In 

2 



A2iV2 



AliVi 



where 



N,N2 



N1 + N2- 1 



(3.35) 
(3.36) 

(3.37) 



Using (2.5), (3.20), (3.26), (3.28), (3.33)-(3.37) we get the relation for the metric (2.1) 
(n = 2) 



Here c 7^ 0, 



for (A^i, N2) = (6, 3), (5, 5), (8, 2) respectively , fi = /i(r, £2) is defined in (3.24) and 

/2 = /2(t, £1, £2, Ci, A^i, N2) = — £l£2^fe(T, £2, Ci), 

iV2 

where functions are defined in (3.12), (3.31) and (3.32). 



(3.38) 
(3.39) 

(3.40) 



4 Exceptional solution for negative curvatures 

Now, we consider the exceptional solution with x = 0. This solution takes place, when Ai, 
A2 < 0, and has the following form 

g = -dts ® di. + E J^^fsg^\ (4-1) 

where D — 1 + Ni+ N2. Here the synchronous-time parametrization is used. Formula (4.1) 
may be obtained by solving eqs. (2.15)-(2.17) and performing the coordinate transformation 
(2.9), (2.10) and appropriate time reparametrization. 
Moreover, it is not difficult to verify that the metric 

g = -dts (^dts + Y. j^zr^t'g^'^^ (4-2) 

defined on the manifold (1.1), where [Mi^g'^''-^) are Einstein spaces satisfying (2.2) (Aj 7^ 0) 
and dimM = £), is a solution to the vacuum Einstein equations (or, equivalently, is Ricci- 
fiat). (See (2.3), (2.4)). 
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Remark 1. We note that the metric 

^•=g(^j'- (4.3) 

defined on the manifold Mi x . . . x M„ satisfies the same relation Rmn[g*] = —{D — 2)g^^ 
as the metric g{H^~^) of {D — l)-dimensional Lobachevsky space. We also note that the 
metric 

= -dt, (g) dt, + tlg{H^-^) (4.4) 

is flat. For D = 3 (4.4) coincides with the well-known Milne solution I^Sf. So, the metric 
(4.2) may be called a "quasi-flat" or "Milne-like" metric. 

Remark 2. It may be proved that, when (Mj, (?(*)) are spaces of constant (non-zero) 
curvature, i = 1, ...,n, and n > 2, the metric (4.2) has a divergent Riemann-tensor squared 
as ts —>■ +0 

m = RMNPQ[g]R''''''%] - +00. (4.5) 

In this case the solution (4.2) is singular. The relation (4.5) follows from the formula 

I[g] = At;\ (4.6) 

where 

^= (D- 2)2^ -2(D -!)(/}- 2). (4.7) 



see 



1^,113]). For spaces of constant negative curvature 



^-2(fl-2)'(|:^-§^)>0, (4.8) 

if n > 2. 

Isotropization. Now, we show that the solution (4.1) is an attractor for solutions (3.38) 
with £1 < 0, £2 < and c > 0, as r — >• 1 (or — > +cxo). Using the relation /2(t) ^ 1 as 
r — > 1, we get the asymptotical formula for the synchronous-time variable tg (c > 0) 

'^"^*.~^§Pl--ir'=-r^l. (4.9) 

Hence 

(/i/2)-^ - -j^tl = as r ^ 1. (4.10) 

Thus, we obtain the relation for the metric 

g = -dts ® dts + h{ts)\Xi\g^^^ + b2{ts)\X2\9^^\ (4.11) 

where 

kits) ~ as ts +00. (4.12) 
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5 Kasner-like behaviour for ts +0 

Let us consider the solutions (3.38) written in the synchronous-time parametrization 

g = -dts ® dts + al{ts)g^^^ + al{ts)g^'^\ (5.1) 

where 

ts = ± dr\hhr''i-' [cj^J + to., (5.2) 

{fi = fi{T'))- As follows from ||13| almost all solutions for the considered model should have 
a Kasner-like behaviour for small t^,, i.e. 

aiiQr^Aie/ ast, ^+0, (5.3) 

where > are constants and the Kasner parameters ctj satisfy the relations 

2 2 

Y,N,a, = Y.N,al = l. (5.4) 

i=l i=l 

Here a suitable choice of the constant tos = ^os('^o) and the sign in (5.2), such that v = 
af^^a^^ — i> +0 corresponds to tg +0, is assumed. We note, that although in [|13| models 
with minisuperspace dimensions n > 3 were considered, the results are also applicable for 
the case n = 2. 



Solving the equations (5.4) we get |Tl 



^ N,±R ^ N2TR , . 

{N, + N2)N,' {N, + N2)N2' ^ ' 

where R = R{Ni, N2) is defined in (2.24). For (A^i, N2) = (6, 3), (5, 5), (8, 2) R is integer 

R = Ni + 2N2 = 12, 15, 12 (5.6) 
correspondingly. We obtain from (5.5) 

^"^'"^)-=G'4)' (14)' ^'-^^ 

«2)- = (-1 1) , (-1 1) , ^) , (5.8) 



9 9/ V 5 5/ V 20 10 

for {Ni,N2) = (6,3), (5,5), (8,2) respectively. Here {ai,a2)± corresponds to "±" in (5.5). 

Now, we return to the metric (3.38) written in the synchronous-time form (5.1). The 
Kasner-like behaviour (5.3) with ctj from (5.7) takes place if r — >• r*, where /2('7"*) = and 
f'li'^*) 7^ 0- Indeed, from (5.2) we get 

ts ~ const|r — r*| 2 as r — > r^, (5.9) 

(r is defined in (3.39)) and using (3.38) we obtain the Kasner-like behaviour (5.3) with the 
parameters 

1 — r — r 

Oil = , 02 = 7^ , (5.10) 

2 — r 2 — r 
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coinciding with those from eq. (5.7) 

The Kasner-hke asymptotes (5.3) with parameters from (5.8) take place, when r —>■ ±oo 
, for non-exceptional values of the constant Ci. 

Ci ^ 0, if 62 = -1 and ^ T^, if 62 = +1, for (iVi, N2) = (6, 3), (8, 2); (5.11) 
Ci ^ ±2 for (ATi, A^a) = (5, 5). (5.12) 

In this case 

/2(r) ~ Eilrl'^ as r ^ ±00, (5.13) 
where B± are constants (integer k is defined in (3.29)) and 

ts ~ const|r|^"^~2 as r ^ ±00. (5-14) 

Here s = {k + 2)r. From (3.38), (5.13) and (5.14) we get the Kasner-like behaviour (5.3) 
with parameters 

«i = |^, «, = ^, (5.15) 

coinciding with (5.8). 

According to the results of |16] the Riemann tensor squared for the solutions with asymp- 
totically Kasner-like behaviour for tg — > +0 is divergent 

I[g] = RMNPQ[g]R^'''''%] - +00 (5.16) 

as ts — »■ -|-0 for fixed {xi,X2) G Mi x M2. If /[fi'*-*-'] > Q for some q, i = 1, 2, then the relation 
(5.16) takes place uniformly on {xi,X2) € Mi x M2 (see Theorem in ||T6[|). 



6 Special solutions with n < 5 curvatures 

The obtained above solutions may be used for generating some special (two-parametric) 
classes of the solutions to vacuum Einstein equations. This may be done using the "curvature- 
splitting" method described below. 

Let us consider a set of k Einstein manifolds (Aii, h^^^) with non-zero curvature, i.e. 

Ric(/i») = (6.1) 

where /Xj 7^ is a real constant, i = l,...,k. Here and below we denote by Ric(/;.) the 
Ricci-tensor corresponding to the metric h. Let 7^ be a real number. Then 

h = j2 ^h^'^ (6.2) 
i=i 

is an Einstein metric, (correctly) defined on the manifold 

M= MiX ...X Mk (6.3) 
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and satisfying 

Ric(/i) = fxh. (6.4) 

Indeed, 

Ric(/i) = Ric(-/i^'^) = E Ric(^^'^) = E /^^^^'^ = (6-5) 

i=l A* i=i i=i 

Remark 3. In (6.2) (like in (2.1)) we identify the metric /i^'^ on M.i with its canonical 
extension to the manifold M. (6.3). It is more correct to write (instead of (6.2)) 

h = jz^plh^^ (6.6) 
1=1 ^ 

where pi : M ^ Mi is the canonical projection. Analogously, Ric(/i^*^) in (6.5) should be 

understood as p* Ric(/i(^)) etc. 

Now, using the suggested trick, we may consider the following Einstein spaces (Mj,gf(*)) 
(i = 1,2) in (3.38): 

gii) = J2 (6.7) 
i = l 

Mi^ Max ...X Min„ (6.8) 
where (M^ , g'-'^^) are Einstein spaces of non-zero curvature 

Ric(^(*^)) = Ai,-^(^^), Kjj^Q, (6.9) 
j — 1, ... ,ni; i — 1,2. Clearly that 



TV.^^AT.. (6.10) 



where 



Nij = dimMi^- > 1, (6.11) 
j = 1, . . .n^; i = 1, 2. It follows from (6.10), (6.11) that 



l<ni< 



(6.12) 



i = 1,2; where [x] denotes the integer part of x. 

Substituting (6.7) into (3.38) we get the following solutions to (vacuum) Einstein equa- 
tions 



9 = cihhr { -l^fi'dT ®dT + Y. heiXi.g^''^ + E^2A2.^(^') ) , (6.13) 

j=i ;=i 



defined on the manifold 

M = X Mil X ... X Mi„j X M21 X . . . X M2n2- (6-14) 
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The solution (6.13), (6.14) describes the evolution oin — ni+n2 Einstein spaces (Mj^, g^'^^^), 
satisfying (6.9)-(6.12). Here {N,,N2) = (6, 3), (5, 5), (8, 2); c 7^ 0, = ±1; r = r(A^i,A^2) 
and fi = fi{T) are defined in subsections 3.1 and 3.2, i = 1, 2. 
The relations (6.13), (6.14) give us 

(1 + 4 + 2 + 3), (1 + 3 + 3 + 3), (1 + 2 + 2 + 2 + 3) (6.15) 

solutions for {Ni,N2) = (6,3); 

(1 + 5 + 3 + 2), (1 + 3 + 2 + 3 + 2) (6.16) 

solutions for (A^i, A^2) = (5, 5) and 

(1 + 6 + 2 + 2), (1 + 5 + 3 + 2), (1 + 4 + 4 + 2), (6.17) 
(1 + 4 + 2 + 2 + 2), (1 + 3 + 3 + 2 + 2), (1 + 2 + 2 + 2 + 2 + 2) 

solutions for {Ni,N2) = (8,2). 

In the last case ((8, 2)) we have two families (of solutions) of the same type as in (6.16). 
(The corresponding solutions from (6.16) and (6.17) seem to be different.) Thus, here we 
obtained some special two-parametrical families of solutions to the problem of cosmological 
evolution of n curvatures for n = 3, 4, 5. 



7 Non-singular solutions 



Here we show that there exist non-singular solutions among the considered ones. Let us 
restrict ourselves to the (1 + 6 + 3)-case, with two negative curvatures ei = S2 = —1 and 
Ci = 0. In this case the functions g{T), 



f2{T) = 3(r^ - 1) 



1 + -rln 
2 



r - 1 



r+ 1 



+ 1 



(7.1) 



and fi{r)f2{r) are holomorphic in O = O U {+00}, where Q = C \ [—1, 1] {Q is the complex 
Riemann sphere with a cut). 



9{r) = - E 



1 



-2fe-l 



2A; + 1 



k=0 
O 00 

k=2 
o 00 

/i(r)/2(r) = - + E/,fcr-^ 
^ k=i 

\t\ > 1. 

Now we introduce a new time variable defined by the relation 



In 



P_ 
Po 



(7.2) 
(7.3) 
(7.4) 



(7.5) 
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where tq > 1, po > 0. We may rewrite (7.5) as 

POTO 



exp[-/(To,T)], 



where 



/(to,t)= rdr' 

J To 



5/1V2 T 



(7.6) 



(7.7) 



The function (7.7) may be analytically continued to the neighbourhood of 00, i.e. in {|t| > 
T} U {00} for some T, and, clearly that 



/(to,t) = /(to,+oo) + ^4(-) . 

k=l 



(7.8) 



We put 

PoTq exp [-/(to, +00)] = 1. (7.9) 

Then p ~ ^ as r — 00 . The function p = p{t) (7.6) is smooth and monotonically decreasing 
on (1, +00) and may be analytically continued to {|r| > T} U {00}, where 



Pir) 



1 + 

k=l 



2k 



(7.10) 



For r — > 1 from (7.5) we get 



p(r) ~ const(T - 1)-Va/io_ 



(7.11) 



Now, we put also 



(7.12) 



where dill is the standard metric on 6-dimensional sphere, normalized by the condition 



Ric {d^fj = bd^l. 
For this special case we get from (3.23), (7.5), (7.12), (7.13) 

g = cFi(p) {-F2(p) [dp ®dp + p'^dnf^ + lAsb^')} , 



where 



i^i(p) - [/l(r(p))/2(r(p))]-^/^ 
F.(p) ^ 



(7.13) 

(7.14) 

(7.15) 
(7.16) 



Here the inverse function r = r(p) is smooth for p e (0, +00) and may be analytically 
continued to the domain {|p| < 5} for some 5 > 0. The same is valid for the functions 
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(7.15), (7.16) (see (2.3) and (2.10)). We have also for \p\ < 5: Fi{p) = Fi{~p), i = 1,2, 
t{—p) = —t{p) and 



-2N -1/2 

Flip) ' 



1 + E Fi,^p"' 

k=l 



(7.17) 



5. 

F,{p)=2 + J2F,,kp'' (7.18) 



k=l 



(see (7.3), (7.4), (7.10)). 

There exist functions $i(w), i = 1,2, defined on {\w\ < 6"^} U (0, +oo), satisfying 

$.(/) = F,(p), (7.19) 

holomorphic in {\w\ < 6"^} and smooth on (0, +cx3). In terms of these functions the solution 
(7.14) may written as 

g = c$i(|xn S^-<^2{\x\^)J2dx' ® dx' + |A2|(7^')| , (7.20) 

where |f p = ELi i.^') ■ Clearly, (7.20) is a smooth metric on the manifold 

X M2. (7.21) 

The solution (7.20), (7.21) is an extension of the special solution defined for r G (0, +00) 
to the semi-interval r G (0, +00]. It may be interpreted in two ways. First, we may say 
that the metric (7.20) describes an extension of the 4-dimensional R x M2 cosmological 
solution to the case of 7-dimensional Euclidean time manifold. On the other hand it may 
be also considered as a spherically-symmetric (0(7)-symmetric) solution with a curved time 
manifold M2 (of negative curvature). 

It seems likely that the procedure considered here may be also applied to the special 
solution with Ci = |, ei = —1, 62 = +1 for A^i = 6, iV2 = 3 and for other special solutions 
in (1 -|- 8 -|- 2)- and (1 + 5 + 5)-models. We may also obtain the "Milne-like" solution, when 
{Mi,g^^^) is the 6-dimensional Lobachevsky space {H^, g{H^)). 



8 Discussion 

In the considered paper we have integrated the vacuum Einstein equations for the model 
describing the evolution of two Einstein spaces Mi and M2 with dimensions {Ni,N2) = 
(6, 3), (5, 5), (8, 2). To our knowledge these are the first non-trivial cosmological solutions 
describing the evolution of more then one Einstein spaces of non-zero curvatures. 

The Kasner-like behaviour near the singularity (for tg — > +0) is investigated. The Kasner 
parameters (5.7), (5.8) are rational for all considered three cases. We may consider the 
following hypothesis: the Abel equation (2.20) may be integrated by methods described in 



20| , pT| (or by its extension) for the dimensions satisfying 



R{Ni, N2) = jNiN2{Ni + N2-l)=meZ. (8.1) 
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(m is integer). In this (and only in this case) the Kasner parameters (5.7), (5.8) are rational. 
(The relation (8.1) is satisfied for (A^i, N2) = (13, 13), (25, 3), (25, 25), (41, 41) etc). 

Here we also received some special solutions to Einstein equations, i.e. "Milne-type" 
solutions for arbitrary n (see (4.2)), and two-paramctric families for 3 < n < 5 (see (6.13)). 
We also obtained non-singular solutions with topology x M2 (sec. 7). (It may be shown 
that the non-singular solutions with topology R^^'^^ x exist also for other dimensions 
(7V,,7V2) = (5,5),(8,2).) 

The considered in this paper solutions to the vacuum Einstein equations are defined on 
the manifold M of dimension D =dimAf = 10, 11. These solutions satisfy the equations of 
motion for D = 10 supergravity of superstring origin and for D = 11, N = 1 supergravity 
respectively. For certain manifolds Mi , M2 it is possible to use the obtained here solutions for 
generation of other classical solutions with non-zero matter fields (dilatonic, Kalb-Ramond 
etc) in D = 10, 11 supergravities. (This may be done, for example, by using the duality 
transformations.) 
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Fig. 1. The graphical representation for the function 2$(r, —S, Ci) (3.12) for S — +1 and 
a) Ci > 0, b) Ci = 0, c) Ci < 0, respectively. 



Fig. 2. The graphical representation for the function 2$(t, — 5, Ci) (3.12) for S — —1 
and a) < Ci < f , b) Ci > f . 

Fig. 3. The trajectories on the plane of scale factors = ai(r),i = 1,2, (corresponding 
to the solution (3.23)) with two negative curvatures (^i = £2 = —1)- Here a) r e (1, +00), 
for Ci > 0; T e (l,T3(Ci)) for Ci < 0; b) r e (t2(Ci), 1). 

Fig. 4. The trajectories on the plane of scale factors aj = ai{T),i = 1,2, for the case 
£1 = +1, £2 = —1, (6- dimensional space has positive curvature and 3-dimensional space 
has negative curvature). Here a) r e (t3(Ci),+oo) for Ci < 0; b) r e (ti(Ci), T2(Ci)) for 
Ci > 0. 

Fig. 5. The trajectories on the plane of scale factors = ai{T),i = 1,2, for the case 
£2 = +1 (3-dimensional space has positive curvature) and a) Si — +1 (6-dimensional space 
has positive curvature), r G (— oo,ri(Ci)) for Ci > |; b) ei = —1 (6-dimensional space has 
negative curvature), r e (— oo,-|-oo) for Ci < |; r e (ti(Ci), -|-oo) for Ci > |. 
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